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Abstract 

The properties and origin of defect modes in the one 
dimensional photonic crystal (PC) with a central defect has 
been studied. Two types of photonic crystals, each having 
single defect, namely, symmetric and asymmetric PCs are 
considered. It is found that a defect mode arise at central 
wavelength for asymmetric PC, whereas two defect modes 
arise in the vicinity of central wavelength for symmetric PC. 
These two defect modes can be fixed to a single central 
defect by increasing the width of the defect layer. But the 
intensities of these defect modes are found to be different 
from the intensities of defect modes for asymmetric case. The 
origin of these defect modes can be explained using 
impedance matching condition. The propagation 
characteristics of the proposed structure are analyzed by 
using the transfer matrix method. 
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Introduction 

Photonic crystals (PCs) which are artificial structures 
with periodically modulated dielectric constants have 
been attracting a great deal of interest among the 
researchers particularly for the study of their 
electromagnetic properties [1-7]. It was observed that 
periodic modulation of the dielectric functions 
significantly modifies the spectral properties of the 
electromagnetic waves. The electromagnetic spectrum 
in such structures is characterized by the presence of 
allowed and forbidden photonic energy bands similar 
to the electronic band structure of periodic potentials. 
For this reason, such a new class of artificial optical 
material with periodic dielectric modulation is known 



as photonic band gap (PBG) material [8]. Fundamental 
optical properties like band structure, reflectance, 
group velocity and rate of spontaneous emission, etc. 
can be controlled effectively by changing the spatial 
distribution of the dielectric function [4, 5]. Photonic 
crystal structure has many interesting applications in 
the field of photonics and optical engineering. 

PCs can guide the flow of light inside them, and they 
have many applications such as localization of the 
light wave [9], inhibition spontaneous emission [10], 
lasers [11-13], waveguides [14], splitters [15], antennas 
[16], optical fibres [17], ultrafast optical switches [18], 
optical circuits [19], tuneable optical filters [20] and 
absolute omni-directional PBGs [21-25]. 

These applications can be realized using pure PCs, but 
doped or defective PCs may be more useful, just as 
semiconductor doped by impurities are more 
important than the pure ones for various applications. 
The idea of doped PCs comes from the consideration 
of the analogy between electromagnetism and solid 
state physics, which lead to the study of band 
structures of periodic materials and further to the 
possibility of the occurrence of localized modes in the 
band gap when a defect is introduced in the lattice. 
These defect-enhanced structures are called doped 
photonic crystals and present some resonant 
transmittance peaks in the band gap corresponding to 
the occurrence of the localized states [26], due to the 
change of the interference behaviour of the incident 
waves. Defect(s) can be introduced into perfect PCs by 
changing the thickness of the layer [27], inserting 
another dielectric into the structure [28], or removing a 
layer from it [29, 30]. 
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PCs with defect(s) have received attention of a huge 
number of researchers, especially 2D and 3D Pes, 
because a great number of applications that can be 
performed using them. It has been also known that a 
point defect in 2D PCs can act as a micro-cavity, a line 
defect in 3D PCs behaves like a waveguide and a 
planer defect in 3D PCs behaves like a perfect mirror 
[31, 32]. Similar to 2D or 3D PCs, the introduction of 
the defect layers in ID PCs can also create localized 
defect modes within the PBGs. Due to the simplicity in 
ID PCs fabrications as compared to 2D and 3D PCs, 
the defect mode can be easily introduced in ID PCs for 
various applications such as in the design of optical 
filters and splitters, in the light emitting diodes and in 
the fabrication of lasers [33-35]. 

In the present communication, the properties of defect 
modes for different combinations of materials and the 
origin of these defect modes have been theoretically 
studied by using transfer matrix method. The origin of 
these defect modes can be explained using impedance 
matching condition. 

Theoretical Analysis 

The schematic representation of one- dimensional 
photonic crystal with defect is represented in Figure 1. 

We consider two structures (one asymmetric and other 
symmetric) in which H and L represents the high and 
low refractive index materials. To compute the 
transmission spectrum, we employ the transfer matrix 
method (TMM) [36]. In this method, the transfer 
matrix for each layer can be written as 



M = F.P.FT 

j j j j 



(1) 



where, j stands for H or L layers and Fj and Pj are 
called the dynamical matrix and the propagation 
matrix respectively. The dynamical matrix is given by 
the following equations 
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FIG. 1 SCHEMATIC DIAGRAM OF 1-D PHOTONIC CRYSTAL (A) 
ASYMMETRIC STRUCTURE (B) SYMMETRIC STRUCTURE 



where the phase is written as 
7.nd , 
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The transfer matrix, for the two PCs can be written as 

( M u M 12 
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\M V M 
for asymmetric PC 
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for symmetric PC 
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The reflection and transmission coefficients in terms of 
the matrix elements, given by Equations (6) and (7), 
can be written as 



and 



_ (M u +q f M l2 )q i - (M 2l +q f M 22 ) 
(M n +q f M 12 )q t +(M 2l +q f M 22 ) 
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(9) 



where qu=ni,fcos9j,f for TE wave and qy= (cos0j,f)/ny for 
TM wave, where the subscripts i and f correspond to 
the quantities respectively in the medium of incidence 
and the medium of emergence. Whereas, the 
reflectance and transmittance of the structure is given 
by 



R = r and y 



<li 



(10) 



In the next section, we will compute the transmission 
spectra of photonic crystal with defect by using 
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Equation (10) and are represented graphically in 
Figure 1. 
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FIG. 2 TRANSMBSDN SPECTRA OF (A) AIR/(HL) 16 /SUBSTRATE 
(B) AIR/(HL) 8 D(HL) 8 /SI02 STRUCTURE 

Results and Discussion 

Let us consider a conventional PC- air/(HL)i6/Substrate 
that has a Photonic band gap in the ultraviolet and 
visible regions. In our study, H is taken to be TiCh 
with refractive index (n) 2.34867 and L is taken to be 
MgF2 with n=1.3855. MgF2 is the good material for the 
PCs because its absorption begins below 115nm and 
the real part of the refractive index is close to the unity 
[37]. We use SiCfe (n=1.46) as the substrate. The 
thickness of each layers is taken according to quarter- 
wavelength condition i.e. nHdH=nLdL=Ao/4, where Aois 
the design wavelength and taken to be 350nm. We 
limit our study to normal incidence only. We have 
plotted the transmittance properties of this structure 
using Equation (10). The transmission spectrum of this 
PC is shown in Figure 2(a). From this figure, it is clear 
that a PBG is found in the wavelength range of 305- 
420nm in the vicinity of Ao. According to the theory of 
dielectric mirrors, the lower and upper band edges of 
the band gap can be calculated using the relation 
given in relation (11) [38, 39]. 

cos 2 f *("» d » + j - k 2 cor f *("» d « - r hM^ \ (11) 

where k=(nH-nL)/(nH+nL) is given by the Fresnel 
formula. By solving Equation (11), we can find the 
lower and upper edge of the PBG which are AL=305nm 
and AH=420nm. If we increase the refractive index 
contrast (nH/m), the width of the PBG increases 
accordingly. If we introduce a defect SiCh (n=1.46) 



layer in the middle of this conventional PC, then a 
defect mode is found at the designing wavelength Aoin 
the PBG as shown in Figure 2(b). Our aim is to study 
this defect mode for various materials namely SiCh, 
AI2O3, Bi4Ge30i2 (BGO), ZnS. We consider asymmetric 
and symmetric PC structures for our study. Also, the 
origin of these defect modes explained using 
impedance matching condition. 
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FIG. 3 TRANSMBSDN SPECTRA OF AIR/(HL) 8 D(HL) 8 /SIQ2 
STRUCTURE (Dd= Ao/4Nd) FOR (A) EMSIO2 (B) EfeALzOs (Q 
D=BGO (D) D=ZNS 
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FIG. 4 PLOT OF T (EQUATION (10)), X AND Y (EQUATION (12)) 
IN THE VICINITY OF THE DEFECT MODES OF FIGURE 3 

Effect of the Defect Material on Defect Mode 

(1) Asymmetric PC structure air/(HL)s/D/(HL) 8 /Si02 

For this case we choose the materials whose refractive 
indices are same as that of H and L layers considered 
for the conventional PC. For the defect layer D, we 
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choose four different materials SiCh (n=1.46), AI2O3 
(n=l .66574), BGO (n=2.13) and ZnS(n=2.58789) [40]. 
The thickness of the defect layer is also chosen so that 
quarter wavelength condition i.e. dD= Ao/4nD is 
satisfied by the defect layer. The transmission spectra 
for these structures, each with a defect, are shown in 
Figure 3(a-d). From these figures, it is clear that for all 
structures with defect of different materials, the defect 
modes are found at the same wavelength Ao. But the 
intensities of these defect modes (transmitted) are 
different for different materials. The intensity is 
maximum for the material whose refractive index is 
closer to the refractive index of low refractive index 
material i.e. MgF2. If we increase the refractive index 
of the defect layer, then the intensity of defect modes 
decreases. 
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FIG. 5 TRANSMBSDN SPECTRA OF AIR/(HL) 8 D(LH) 8 /SI02 
STRUCTURE (Dd= Ao/4Nd) FOR (A) D=SICfe (B) ENAL2O3 (Q 
D=BGO (D) D=ZNS 
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FIG. 6 TRANSMBSDN SPECTRA OF AIR/(HL) 8 D(HL) 8 /SI02 
STRUCTURE (Dd= Ao/2Nd) FOR(A) D=SI02 (B) ENAL2O3 (Q 
D=BGO (D) D=ZNS 

The origin of these defect modes can be explained 
using impedance matching condition [38]. The 



normalized effective surface impedance, z at the 
incident plane boundary is given by 



z = X+iY-- 



Z _l + r 



(12) 



where Zo=377Q, is the intrinsic impedance of the free 
space, Z is the effective surface impedance at the 
incident plane boundary and r is the reflection 
coefficient defined by Equation (8). The real part X, 
imaginary part Y (Equation (12)) and transmittance, T 
are plotted in Figure 3(a-d) in the vicinity of defect 
modes. From these Figures, it is seen that the defect 
mode is observed only when the condition of 
impedance matching is satisfied. In these Figures, we 
have found an anomalous behaviour in real part X of 
the effective surface impedance Z. The intensity of 
defect modes is found be maximum corresponding to 
the strong dispersion of X and it is, however, be 
minimum corresponding to weak dispersion of X. 

(2) Symmetric PC structure air/(HL)8/D/(LH)s/Si02 

For this case, we choose the same materials as taken in 
case (1). The transmission spectra for symmetric 
structure for different materials are shown in Figures 
5(a-d). From these Figures, it is clear that in this case, 
two defect modes are found near the band edges of the 
PBG. These defect modes disappear if we increase the 
refractive index of the defect layer (Figure 5(d)). 
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FIG. 7 PLOT OF T (EQUATION (10)), X AND Y (EQUATION (12)) 
IN THE VICINITY OF THE DEFECT MODES OF FIGURE 6 

On the other hand these defect modes can be fixed to a 
single defect mode at the design wavelength (Ao) by 
increasing the width of the defect layer twice its 
previous value i.e. dD= Ao/2nD. The transmission 
spectra for this case are shown in Figure 6(a-d). From 
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these Figures it is clear that the intensity of each of 
these defect modes is more for those defect materials 
whose refractive index is close to the high refractive 
index material (H). For this case the plot of X, Y and T 
is shown in Figure 7(a-d). From these Figures it is clear 
that there is a strong dispersion corresponding to a 
high intensity of the defect mode. 




275 350 425 
Wavelength (run) 



350 425 
Wavelength (nm) 






< 


i. 








1 






! 










I ji 
•. -i 




i < 






; -i 










■r .i 








i 


! ' ■' 








i : 






d) 









275 350 425 
Wavelength (nm) 



:::: 



275 350 425 
Wavelength (nm) 



FIG. 8 TRANSMBSDN SPECTRA OF AIR/(HL) S D(HL) 8 /SI02 
STRUCTURE (Dd= Ao/4Nd) FOR (A) D=SIQ> (B) D=AL2Qi (C) 
D=BGO (D) D=ZNS AT DIFFERENT ANGLES OF INCIDENCE CP 
(BLACK SOLID LFNE), 30 a (BLUE DOTTED LINE) AND 60 s (RED 
DASHED LINE) 




1 



ML 



(a) 



275 350 425 
Wavelength (nm) 



KZ 




275 350 425 
Wavelength (nm) 



uhJitlf li 








ill 










ft* 'M 

5 .'■«! -:» 


,i" : 
i'i. 






'■!<. 




J» :;| 


i ■■"'i 




i :■: 






* ■■■ 1 










(c) 




1 ■■*■.! 



13 

Ed5 





| J 

N 


fc'.l I'J I 


Mm 








H 1 

VI 

' l' 


i : ■ 





w 



350 425 
avelength (nm) 



275 350 425 
Wavelength (nm) 



FIG. 9: TRANSMBSDN SPECTRA OF AIR/(HL) 8 D(LH) 8 /Sia 
STRUCTURE (Dd= Ao/4Nd) FOR (A) D=SI02 (B) ENAL2O3 (Q 
D=BGO (D) D=ZNS AT DIFFERENT ANGLES OF INCIDENCE s 
(BLACK SOLID LINE), 30 a (BLUE DOTTED LINE) AND 60 s (RED 
DASHED LINE) 




275 350 425 
Wavelength (nm) 



275 350 425 
Wavelength (nm) 



FIG. 10: TRANSMBSDN SPECTRA OF AIR/(HL) 8 D(LH) 8 /SI02 
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Effect of Angle of Incidence on Defect Mode 

(1) Asymmetric PC structure air/(HL)8/D/(HL) 8 /Si0 2 

In case of the asymmetric structure, the transmission 
spectra for different defect layers SiCh, AI2O3, BGO and 
ZnS with defect thickness dD= Ao/4nD, at different 
angles (0 2 , 30° and 60 B ) of incidence are shown in 
Figures 8. From this Figure, it is clear that the defect 
modes shift towards the lower edge of the band gap as 
we increase the angle of incidence. Also the intensity 
of these defect modes decreases with increase the 
angle of incidence. It is also clear that the intensity of 
these defect modes decreases with fewer amounts for 
those defect materials whose refractive index is close 
to the (nH~nL,)/2. 

(2) Symmetric PC structure air/(HL) 8 /D/(HL)s/Si02 

In case of the symmetric structure, the transmission 
spectra for different defect layers SiCh, AI2O3, BGO and 
ZnS with defect thickness dD= Ao/4nD, at different 
angles (0 s , 30 9 and 60 s ) of incidence are shown in 
Figure 9. From these Figures, it is clear that in this case, 
two defect modes are found near the band edges of the 
PBG at normal incidence. As we increase the angle of 
incidence the left defect mode disappear whereas the 
right defect mode shifts towards the left. Both defect 
modes disappear with the increment of the refractive 
index of the defect mode. Also the intensity of these 
defect modes decreases with increase the angle of 
incidence. These defect modes disappear if we 
increase the refractive index of the defect layer. 
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On the other hand these defect modes can be fixed to a 
single defect mode at the design wavelength (Ao) by 
increasing the width of the defect layer twice its 
previous value i.e. dD= Ao/2nD. For this case of the 
symmetric structure, the transmission spectra for 
different defect layers SiCh, AI2O3, BGO and ZnS at 
different angles of incidence (O 2 , 30 s and 60 2 ) are 
shown in Figures 10. From this Figure, it is clear that 
the defect modes shift towards the lower edge of the 
band gap as we increase the angle of incidence. Also 
the intensity of these defect modes decreases with 
increase the angle of incidence. In this case, the 
intensity of the defect modes decreases in the different 
way from the asymmetric case. 

Conclusions 

Based on the Transfer Matrix Method, the properties 
and origin of defect modes in one dimensional 
photonic crystal have been theoretically studied. It is 
found that we can achieve the desired intensities of 
defect modes using symmetric and asymmetric 
structures and the thickness of the defect layer 
according to our requirement. These defective PCs can 
be used as single and multi channel filters. Specially, 
the symmetric defective PC can be used also in 
DWDM applications. 
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